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Abstract 

The properties of open charm mesons, D, D, D s and D s in nuclear matter at finite temperature are 
studied within a self-consistent coupled-channel approach. The interaction of the low lying pseudoscalar 
mesons with the ground state baryons in the charm sector is derived from a i-channel vector-exchange 
model. The in-medium scattering amplitudes are obtained by solving the Lippmann-Schwinger equation at 
finite temperature including Pauli blocking effects, as well as D, D, D s and D s self-energies taking their 
mutual influence into account. We find that the in-medium properties of the D meson are affected by the 
-D s -meson self-energy through the intermediate D S Y loops coupled to DN states. Similarly, dressing the D 
meson in the DY loops has an influence over the properties of the D s meson. 

PACS numbers: 14.20.Lq, 14.40. Lb, 21.65.+f, 12.38.Lg 
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I. INTRODUCTION 



Over the past years the properties of charmed hadrons have received a lot of attention in 
connection with experiments in lepton colliders (CLEO, Belle, BaBar) and hadron facilities (CDF 
ab, PHENIX, STAR at RHIC, and the forthcoming PANDA and CBM experiments at 



at Fermi 

FAIR) ll422l] . The discovery of new resonances with charm content has sparked the interest of not 
only many experimental but also theoretical research groups in order to find plausible explanations 
for the nature of the newly found states. 

The study and characterization of resonances has been a very active topic of research in hadrons 
physics during the last decades. The goal is to establish whether some resonances have the genuine 

I24I or, alternatively, qualify better as hadron 



qq or qqq structure predicted by the quark model [2 
molecules generated dynamically. A series of pioneer works 25l430l|. based on a t-channel vector- 
meson exchange (TVME) force, and more recent approaches in terms of chiral Lagrangians 



42 



44 



491 ] have proven to be very promising and successful in describing a wealth of S- wave baryonic 



resonances in the light SU(3) sector within the molecular picture using coupled-channel dynamics 
with effective hadronic degrees of freedom. In the modern language of chiral Lagrangians those 
states emerge from the scattering of the 0~ octet Goldstone bosons off baryons of the nucleon 
l/2 + octet. Moreover, P— and D— wave baryonic resonances have been generated based on 



the SU(3 
decuplet 



leading-order chiral Lagrangian by means of meson scatte ring off baryons of the 3/2 + 



50H53I] and incorporating vector-meson degrees of freedom [54-59]. Molecular states of 



two pseudoscalar mesons and one baryon 



6014641 have also been studied to interpret low-lying 



1/2+ states. All these results sustain the so-called hadrogenesis conjecture, formulated by Lutz 
and Kolomeitsev a few years ago, according to which resonances not belonging to the large -/V c 
ground state of QCD are generated by coupled-channel dynamics Q 5 Q]- 

Recently, the charm degree of freedom has been incorporated in those unitarized coupled-channel 
approaches to describe open- and hidden-charm mesons 



in the meson-baryon sector 



71 



6JH70J]. Similar methods have been used 



751 ]. partially motivated by the parallelism between the A(1405) 



and the A c (2595). The meson-baryon interaction in the charm sector is constructed using the t- 
channel exchange of vector mesons between pseudoscalar mesons and baryons and performing the 
zero-range approximation while preserving chiral symmetry for light mesons 73J] . The extension to 



D-wave J = 3/2 resonances was first attempted by extending the basis to include the J = 3/2 + 



baryons [76]. In order to be consistent with the spin-flavor heavy-quark symmetry (HQS) 



the vector mesons were incorporated later within a SU(8)-inspired model 



80 



77 



Z9|, 



811 ] . similar to the 



2 



SU(6) one developed in the light sector |55l. l56l|. An alternative approach based on the local hidden 
gauge formalism has recently become available [sj j. On the other hand, there have been also 
attempts to construct the DN and DN interaction by incorporating the charm degree of freedom 



83|485]. 



in the SU(3) meson-exchange model of the Jiilich group 

Nuclear medium modifications have been lately incorporated as a second step. The aim is to 
further investigate on the nature of resonant states, such as A c (2595), but also to test the dynamics 
of charmed hadrons with nucleons and nuclei. The properties of open-charm mesons in nuclear 
matter can influence the charmonium production in hot dense matter, which might indicate the 
formation of the quark-gluon plasma phase of QCD at high densit y a nd temperature 
exciting scenario is the possible formation of D-mesic nuclei 
states like J in nuclei 



87|, 



861 ] . Another 



881 ] and of exotic nuclear bound 



89-91]. From the experimental side, the physics program of the CBM 
experiment as well as part of the PANDA collaboration at FAIR [92] will be devoted to the 
properties of open and hidden charm in dense matter. In particular, the physics goal is to extend 
to the heavy-quark sector the GSI program for in-medium modifications of hadron properties in 
the light sector, and to provide insight into the charm-nucleus interaction. 

Works based on mean-field approaches provided important shifts for the D and D open-charm 
meson masses [93j-[96|], wh ich a lters the formation of charmonium [97]. Some of those models have 



been recently revised 



1001 ] . A different perspective is offered by models that, working within 
coupled-channel unitarized schemes, go beyond mean field and p rovide the spectral feat ures of th e 
charm mesons in symmetric nuclear matter at zero 



1011 ] and finite temperature 



102, 



1031 ]. 



Lately, this meson-baryon basis has been extended to incorporate HQS. In this way, not only 



Z}-meson but also D*-meson features have been studied [1041 ] . 



A common feature of the previous models is the use of an interaction kernel in the zero-range 
approximation (t — > 0). This is justified for diagonal amplitudes close to threshold and for non- 
diagonal transition amplitudes where the masses of mesons and of baryons in the initial and final 
meson-baryon states differ moderately. However, the charm-exchange processes, for which the 
difference in masses between the external mesons is comparable with the mass of the charmed 
vector meson being exchanged, point towards the breakdown of the zero-range approximation. 
Charmed baryon resonances ha ve be en studied using the full t-dependence of the i-channel vector- 



exchange driving term in Ref. 



1051 ]. Compare d to the previous TVME local models, where the 



t — > limit was implemented, the work of Ref. 



1051 ] obtained the same amount of resonances but 



located in general at somewhat higher energies and having larger widths. Some of these resonances 
could clearly be identifiable with experimentally seen states, such as A c (2595), S c (2800), H c (2790) 
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and H c (2980). 

In the present work, we study the behavior of the dynamically-generated baryonic resonaces 
in hot dense matter, as well as the spectral features of the open charm mesons (D, D, D s and 
D s ), within a self-consistent coupled -cha nnel approach that considers the full t-dependent TVME 



interaction kernel employed in Ref. [105(. We pay a particular attention to the influence that the 
dressed mesons exert on each other. We find that the simultaneous dressing of the charm mesons 
(D,D S ) in the C = 1 sector, or the anticharm mesons (D,D S ) in the C = —1 one, affects their 
in-medium properties in a non-negligible way. 

The article is organized as follows. In Sec. HTJ we present the formalism. We first revise the model 
adopted for the free space amplitudes and, next, we describe the modifications that incorporate 
the medium effects. Our results for the medium modified resonances and for the spectral functions 
of the open-charm mesons at various densities and temperatures are shown in Sec. III. A summary 
of our conclusions is presented in Sec. IV. 

II. FORMALISM 

In this section, we will first review briefly the coupled-channel approach employed in our previous 
work 



1051 ]. where we studied open-charm baryon resonances dynamically generated from the free- 
space interaction of the low-lying pseudoscalar mesons with the ground-state baryons using a 
t-channel vector-exchange driving force. After that, we will introduce the main sources of medium 
effects and we will implement them in our coupled-channel formalism. 

Since the properties of the D, D, D s and D s mesons in a hot and dense environment will be 
determined, respectively, from the DN, DN, D S N and D S N amplitudes, we list in Table [J the 
corresponding set of coupled channels in each of the related isospin (/), strangeness (S) and charm 
(C) sectors. 

A. Free-space coupled-channel approach 

The free-space amplitudes, T, which describe the scattering of the pseudoscalar meson fields 
off the ground-state baryon fields can be obtained by solving the well-known Lippmann-Schwinger 
equation, which schematically reads 

T = V + VJT . (1) 
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(I,S,C) 


Channels 


(1,-1,-1) 


D S N(2908) L»A(2985) L>£(3062) 


(0,0,-1) 


Z)A(2806) 


(1,0,-1) 


LW(2806) 


(0,0,1) 


ttS c (2591) DA(2806) 7yA c (2832) £TS C (2963) KE' C {3070) 
L> S A(3085) r/A c (3243) ?7 C A C (5265) L>S CC (5307) 


(1,0,1) 


ttA c (2424) ttE c (2591) L>A(2806) A'S C (2963) ryE c (2999) 
ATS' C (3070) L> s £(3162) r/E c (3410) £>S CC (5307) 7y c E c (5432) 


(1,1,1) 


ifA c (2779) D S A(2908) if£ c (2946) 



TABLE I: Coupled-channel meson-baryon states involved in DN, DN, D S N, or D S N interactions. The 
energy threshold of each state is given in brackets. 



The loop function J is the product of the meson and baryon single-particle propagators, and the 
scattering kernel V describes the interaction between the pseudoscalar mesons and the ground-state 
baryons. Following the original work of Hofmann and Lutz [^J], we identify a t-channel exchange 
of vector mesons as the driving force for the S-wave scattering between pseudoscalar mesons in 
16-plet and baryons in 20-plet representations. The scattering kernel takes the form (see [3] for 
details) 



v}j' S ' C) (Pi,Qi,Pj,qj) 



Ve[i6] 



HIT 



1 



t — m 



V 



-{<n+qj) u u(pi) 

(2) 



where the sum runs over all vector mesons of the SU(4) 16-plet, (p, K* , K* , u, <fi, D* , D*, D* , 
D*, J/*S>), my is the mass of the exchanged vector meson, g is the universal vector meson coupling 
constant, Pi,qi,Pj and qj are the four momenta of the incoming and outgoing baryon and meson, 

a s o 

and the coefficients C\j'y denote the strength of the interaction in the different (I, S, C) sectors, 
and me son- baryon channels The value of g = 6.6 that reproduces the decay width of the p 



meson 



1061 ] has been considered in this work. The S-wave projection of the scattering kernel is 
easily obtained, and in the center-of-mass (cm.) frame it takes the analytical form 

,2 



V: 



(I,S,C), 



ij,l=0 



ve[i6] 



2/3 ab- (3a, ( a + b 

T + ~ b^ ln [—b 



(3) 



5 



with a, b, a and f3 being 



a = m} + m 2 j -2ui{\ki\)uj{\kj\) - m v 



3 

1 1 



1 1 

m^ — m- 



a = fiiflfcl) + - Mi - Mj - J 2 (Qjdfcj-I) - Oid^l) + Mj — Mj) 

% 

(Ei(\ki\) + Mi)(Ej(\kj\) + Mj) { * U lU Al Jl ' 

11 \ 

— '-^(njdkjD-nidhD + Mj-Mi)) , (4) 

m v j 

where ki,kj are the initial and final relative momenta, mi,mj, Mi, Mj are the masses of the in- 
coming and outgoing mesons and baryons, and Ui(\ki\),ujj(\kj\), Ei(\ki\), Ej(\kj\) their correspond- 
ing energies, which have been taken to be their on-shell values. The factor N = [(E(\ki\) + 
Mi)(E(\kj\) + Mj) / {AMiMj)] 1 / 2 comes from the normalization of the Dirac spinors. We have de- 
fined 0(|A;|) = w(|fc|) + E(\k\). We note that the zero-range approximation (i.e., t — > 0) of the 
S*-wave scattering kernel is obtained by expanding the logarithm of Eq. (J3]) in the limit b/a — > up 
to the linear term in b/a and setting a = —m v . The interested reader is referred to our previous 



work of Ref. |l05l ] for a detailed analysis of the validity of the zero-range approximation. 

In Eqs. (|2|) and (|3|), we have assumed infinitely (zero- width) exchanged vector mesons, because 
the value of t is never larger than the square of the minimum energy required for the meson to 
decay. In other words, since the mesons being exchanged in this problem are largely off shell, they 
will be treated as stable particles. 

Once the scattering kernel has been constructed, one can finally write the S-wave projection of 
the Lippmann-Schwinger equation, 

4jf^\k h kj, y/8) = V^\k h kj) + £ | ^F(|£|)<£ C) (4 k) MVS, k)T^\k, kj, yfs) 

(5) 

where \fs is the total energy in the cm. frame. The loop function J explicitly reads 

J {I ' S ' C) (J~s k) = — - (6) 

k ,J 2E k {\k\)u k {\k\)^s-E k {\k\)-u Jk (\k\) + ir } ' 

and is a dipole-type form factor, 

F(\k\) = ( A \ ] , (7) 
\A 2 + \k\ 2 J 

that has been introduced to regularize the integral. This form is typically adopted in studies of 
hadron-hadron interactions within the scheme of Lippmann-Schwinger-type equations in the light 



flavour sector 1071 ]. The value of the cut-off A is a free parameter of our model. Given the limited 
amount of data for charmed baryon resonances, and in order to simplify the analysis, the cut-off A 
is adjusted to 903 MeV/c in order to reproduce the position of the well-known J p = l/2~ A c (2595) 
having (/, S, C) = (0, 0, 1), and the same value is used for the other sectors explored in this work. 
In Table [III we summarize the position, width and most important couplings of the dynamically 
generated states appearing in the various (I, S, C) sectors listed in Tabled Note that there are no 
resonances in the singled-channel (I, S, C) = (0, 0, —1) and (1,0, —1) sectors of the DN interaction. 
In the other C = — 1 case, having (/ = \,S = -1), we find a pole just below the P*N threshold. 



1051 ] . we briefly 



The remaining cases have C = 1 and, although they were deeply analyzed in Ref. 
comment here a few essential characteristics that will be useful for our discussion of the in-medium 
results in the next section. In the (/ = 0, S = 0) sector, apart from the A c (2595) resonance to 
which we fit the model, there is another very narrow one at 2805 MeV, just below the threshold 
for DN states but coupling very little to them. We also predict two narrow resonances in the 
(J = 1, S = 0) sector at 2551 and 2804 MeV, right below the thresholds of the channels to which 
they couple more strongly, namely 7rS c and DN, respectively. In the (J = = 1) case, we 
predict a cusp-like structure placed at the threshold of KY, C , the channel that shows the largest 
coupling to this state. 



(I,S,C) 


M R [MeV] 


T [MeV] 


Main decay modes \g\ {channel) 


(1.-1,-1) 


2906 (pole) 





l.3(D s N), l.l(DA), 1.9(Z)£) 


(0,0,1) 


2595 
2805 


0.5 
0.01 


0.31(tt£ c )*, ll(DN), 6.0(£> s A),2.0(t? c A c ) 
0.04(^£ c )*, 0.27(DV), 2.2(nA c ), 4.3(KE C ), 0.21(I> fl A) 


(1,0,1) 


2551 
2804 


0.16 

5 


0.05(ttA c )*, 3.7(tt£ c ), l.l(DJV), 2.1(KE' C ) 
0.27(vrA c )*, 0.14(7rE c )*, 2.1(DJV), 1.8(D a S) 


(1,1,1) 


2946 (cusp) 


0.93 


0.002(A:A C )*, 0.03(D S N)*, 0.07(ifE c ) 



TABLE II: Dynamically generated baryon resonances with open charm in various (I, S, C) sectors for a 
cut-off momentum A = 903 MeV/c. The table shows the position (Mr) and width (T) of the resonance, 
together with the couplings to the most important meson-baryon channels, as well as the couplings to the 
channels in which it can decay (marked with an asterisk). 
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B. Medium effects 



The are two main sources of medium effects to consider: one is a consequence of the Pauli 
exclusion principle, that prevents the scattering of two nucleons into states which are already 
occupied. The other is related to the fact that the properties of all mesons and baryons are 
modified in the medium due to their interactions with the Fermi sea of nucleons. Pauli blocking 
and finite temperature effects can be incorporated in the coupled-channel equations by simply 
replacing the free nucleon propagator by the in-medium one, 



G N {po,p,p, T) 



En{\p\) 



l-n N (p,p,T) n N (p,p,T) 



(8) 



po - E N (\p\) + ie po - E N (\p |) - ie_ 
where (po,p) is the total four- momentum of the nucleon in the nuclear matter rest frame, 
njsr(p,p,T) is the usual Fermi-Dirac distribution function, and En{\p\) is the on-shell energy 
of the nucleon. 

The nuclear medium effects on the mesons can be incorporated by including their corresponding 
self-energies, H m (qo,q, p,T), in the meson propagator 

D m (q ,q,p,T) = * , (9) 

% 1 ,H m 11 m WO) y> Hi 1 ) 

being (qo,q) the four- momentum of the meson. This is done in practice through the corresponding 
Lehmann representation of the meson propagator 

D m {qo,q,p,T = / — — dw- / ■ —du , (10) 

where S m ^(u!, q, p,T) is the spectral function of the meson m(fh): 

a i ■+ rps 1 r / n / -» rji\ \ 1 lm(U m (q ,q,p,T)) 

Sm(qo,q,P,T) = Im (D m (q ,q,p,T)) = — ^ — = —— — —k . (11) 

vr vr |qg - q l - m l m - Ll m (q ,q, p,T)\' s 

We note here that in this work only the D,D, D s and D s mesons have been dressed by self-energy 



75|, 



103 



104 ]. 



insertions. Mesons 7r, K, r/, rf and r] c have not been dressed, as done e.g., in Refs. 
The reason is that the states containing these mesons couple weakly to the DN and D S N ones and, 
therefore, it is expected that approximating the vr, K, rj, rj' spectral functions by the free-space ones, 
i.e., delta functions, will not influence much the in-medium properties of the D and D s mesons. We 
emphasize, however, that the present work addresses for the first time the simultaneous dressing of 
the D and D s mesons in the charm C = 1 sector, and that of the D and D s mesons in the charm 
C = — 1 one. 

The loop function for the free case given by Eq. ([6]) must now be replaced by the one including 
the medium and temperature effects on the baryon and meson propagators, as given by Eqs. flSJ) 



and (|10p . Using the Imaginary Time (or Matsubara) Formalism 1081 ] we obtain: 



T {I,S,C) (T} t\ Mk I f°° rl Q ( Tit ^ ~ n(xP + k, p,T) + f(co,T) 
J k { p 0i P ik,p,T) = — — — — =,— [ / du S m (uj,yP - k,p,T) — — =r- 

' " " ' ' ' Pq — oj — E k (\xP + k |) + ie 




du> Srn(uJ,yP - k,p,T)- — =- (J2) 

P + u - E k (\xP + k\) -tej 

where P = 9o+-Efc(|p |), P = Q+P and k = yp-xq, with x = M k /(m k +M k ) and y = m k / (m k +M k ), 
are the total energy, total momentum, and relative momentum of the meson-baryon pair in the 
nuclear matter rest frame, n is the Fermi distribution of the baryon and / is the Bose enhancement 
factor of the meson. In practice, given the nuclear densities and temperatures explored in the 
present work, we can set / = for all mesons and n = for all baryons except for nucleons. One 
might argue t hat t he Bose enhancement factor for the pions should not be ignored. However, as 
tested in Ref. |loj |. the DN amplitudes are insensitive to this factor due to the reduced coupling 
to 7rS c states resulting from the heavy mass of the meson exchanged in the transition potential. 

The in-medium scattering amplitudes T are obtained by directly solving the coupled-channel 
Eq. (|5|) with the medium modified loop function Jm ' (Po: P, k, p, T). The in-medium self-energies 
for the D,D,D S and D s mesons are then obtained by integrating the in- medium diagonal scattering 
amplitudes over the nucleon Fermi sea as 

U D{D) (q ,q,p,T) = J ^n(p,p,T) [r^ ( = ) yPo,Ap,T) + 3r^ ( = 1 ) ^(P ,P, /0 ,r)] , (13) 

for D and D, and as 

U Ds(Bs) (q ,q,p,T) = 4 j ^ n (p,p,T) T^=^(J%, P, p,T) , (14) 

for D s and D s . 

Finally, we note that the self-energies Tl m (m = D, D, D s , D s ) must be determined in a self- 
consistent way since they are obtained from the in-medium scattering amplitudes Tdn, Tq n , Tjj 3 n 
and Tjj aN , which contain the loop functions Jjw > ^D^Y (DN case); J^'^' ^ (DN case); Jd'Jn 
(D S N case); and J^'^ , J^y' C ^ (D S N case), that are themselves functions of the self-energies 



III. RESULTS AND DISCUSSION 

We will start discussing our results for the C = 1 mesons, D and D s . First of all, we note 
that their in-medium properties will be influenced by the charm C = 1 baryonic resonances that 
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1=0, S=0, C=l 



S=0, C=l 



nr 



> 


-50 
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P=0,I 


100- 




z 
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150- 


s 





3 vr 
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Dressing D 
Dressing D, D s 



2500 2600 2700 2700 



P [MeV] 



2800 



2900 



3000 



P [MeV] 



FIG. 1: Imaginary part of the 1 = (left panel) and 1=1 (right panel) DN — > DN scattering amplitudes 
in nuclear matter at normal saturation density po and zero temperature, as function of the total energy Po, 
for a total momentum P = and various approximations. 



couple significatively to DN and D S N. From the results of our previous work [1051 ] . summarized in 
Table HH we find the well known A c (2595), coupling very strongly to DN states in the (/ = 0, S = 0) 
sector, and two other resonances, E c (2551) and £ c (2804), coupling also significatively to DN in 
the sector (J = 1, S = 0). The cusp-like structure found in the (J = 1/2, S = 1) sector shows also 
a sizable coupling to D S N states. 

In Fig. [TJ we show the imaginary part of the DN amplitude at normal nuclear matter saturation 
density, po = 0.17 fm -1 , and zero temperature, as a function of the center-of-mass energy Pq, 
covering an energy range that includes the most relevant resonance in each isospin sector, 1 = 
(left panel) and 1 = 1 (right panel), for various approximations. The amplitude in free space (p = 0) 
is also shown (dash-dotted lines), as a reference. When only Pauli blocking effects are included 
(dashed line) we observe that the A c (2595) and E c (2804) states appear displaced to higher energies, 
by about 60 and 50 MeV, respectively. This repulsive effect is well known, and has to do with the 
loss of phase space associated to the fact that the nucleons are forced to occupy empty states that 
are located at momentum states above the Fermi momentum. However, when the dressing of the 
D meson is incorporated self-consistently (dotted line) A c (2595) resonance moves to substantially 
lower energies and the S c (2804) dilutes. This is naturally explained in terms of the D-meson 
strength distribution which, as we will see, shows a quasiparticle peak at a lower energy than 
in free space and a pronounced peak at even lower energies related to A c (2595)A r ~ 1 excitations. 
The reduced in-medium DN threshold opens decay channels for the £ c (2804) which, therefore, 
broadens considerably. As for the A c (2595), its position makes it very sensitive to the low energy 
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strength of the D spectral function and, together with the larger coupling to DN states, explains 
why the resonance acquires such a large amount of attraction. 

In Ref. where the TVME in the t — > limit is employed, a similar behavior is observed 
for the A c (2595). The repulsive shift with respect to the free space position due to Pauli blocking 
effects is compensated by the attractive self-consistent dressing of the D meson. However, the shift 
is smaller in Ref. [3], as it can be seen from Fig. 5 of this reference (model B). Moreover, this 
TVME (t —> 0) model also generates a E c resonance, which lies around 2795 MeV for model B. This 
structure melts down as the dressing of D mesons is incorporated, because of the opening of new 
absortion channels, and stays close to the position wi th on ly Pauli blocking effects, in agreement 
with our present results. In the TVME model of Ref. 101[ both resonances are generated but the 
S c one is localized at a much smaller energy, around 2620 MeV. In that work the self-consistent 
dressing of D mesons results in an attractive shift for both A c and S c dynamically generated 



resonances, in contrast to our results. Within the SU(8) Weinberg- Tomozawa model of Ref. 80], 
where heavy-quark spin symmetry is implemented, the S c state lies around 2900 MeV and has a 
different spin, J = 3/2. In this scheme, the S c resonance couples strongly to the D*N channel 



104]. 



instead of DN and behaves similarly in matter as the A c (2595) 

It is clear from Table HI1 that the strong coupling of the A c (2595) resonance to the D S A channel 
and that of the S c (2804) to D S Y, states, makes it advisable to consider also the medium modi- 
fications of the D s meson in the intermediate D S Y loops of the DN amplitude. This is a quite 
arduous task that, up to our knowledge, has been attempted for the first time in the present work. 
Our results, represented by the solid lines in Fig. [U clearly show the importance of such dressing, 
making the A c (2595) to appear 35 MeV higher in energy with respect to the case of considering 
free D s mesons in the D S Y loops. The changes on the S c (2804) resonance are more moderate. 

The real and imaginary parts of the D and D s self-energies and spectral functions at normal 
nuclear matter saturation density and zero temperature are shown in Fig. [2j as functions of the 
meson energy, qo, and for a meson momentum q = MeV/c. The approximations displayed include: 
Pauli blocking effects; the additional self-consistent dressing of the given meson; and, in the case 
of the D meson, the additional dressing of the D s meson in the D S Y intermediate states coupling 
to DN. 

The features discussed for the DN amplitude in Fig. [1] are also reflected in the imaginary 
part of the D-meson self-energy displayed in the upper left panel of Fig. [2j The middle panel 
shows the corresponding real part of the self-energy, HeH(qo,q = 0), together with the function 

— m? D (thin solid line), such that the crossing points of these two functions signal the appearance 
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FIG. 2: Imaginary (upper panels) and real (middle panels) parts of the D (left panel) and D s (right panel) 
meson self-energies and spectral functions (lower panels), as functions of the meson energy qo, at normal 
nuclear matter saturation density, for q = MeV/c and different approximations. The mass of the D meson 
and the function q 2 , — m 2 D (left panel) and the mass of the D s meson and the function q^ — m 2 D (right panel) 



are shown for comparison. 



of pronounced maxima in the spectral function, as long as the imaginary part does not show a 
pronounced minimum there. Actually, the distribution of the D-meson strength shown in the 
lower-left panel is very rich. All the approximations give a quasi-particle peak located around 35 
MeV below the free Z)-meson mass. In addition, each resonance leaving a signature in the self- 
energy produces a resonant-hole excitation peak in the spectral function, located at a somewhat 
different value of energy due to the complex structure of the self-energy. The common behavior is 
that the resonance-hole modes in the spectral function get displaced such that they move further 
away from the quasiparticle peak. In the case of Pauli blocking, we can clearly distinguish three 
of such modes, associated to £ c (2551)iV -1 , A c (2595)iV~ 1 and S C (2804)7V"- 1 excitations. When 
meson dressing is incorporated, only the A c (2595)A r ~ 1 excitation mode is clearly visible. The 
S c (2804)iV- 1 mode merges with the quasi-particle peak, and the S c (2551)A r one is no longer 
visible in the spectral func tion as compared to th e A n (2595)A^ -1 mode. A similar behavior has been 



observed in Refs. 



75 



1031 ]. In contrast, in Ref. [101|], the S c (2804)A f - 1 mode 
lower energy and mixes with the A C (2595)A^ _1 one, while the quasiparticle peak of the D meson 
experiences a repulsive shift of 32 MeV. It is also worth mentioning that in the SU(8)-inspired 
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model of Ref. 1041 ] the quasiparticle peak appears at slightly lower energies than the free mass but 
the D-meson spectral function shows a completely different shape due to the different resonant-hole 
composition of the D-meson self-energy. 

The imaginary part of the D s self-energy, displayed in the upper right panel of Fig. [2] shows 
only a small enhancement at around 2 GeV. This is a refle ction of the enhanced cusp found in the 



(I = 1/2, S = 1,C = 1) amplitude at the KT, C threshold 1051 ] . This structure generates a small 



but non-negligible amount of strength in the D s spectral function to the ri ght of the quasi-particle 



1011 ] , where a resonance 



peak, which barely moves from its free location. This is in contrast to Ref. 
is generated dynamically 75 MeV below the D S N threshold, and the corresponding resonance-hole 
state in the spectral function appears on the left-hand side of the quasiparticle peak. 

In spite of the featureless aspect of the D s spectral function in our model, this relocation of 
strength from the quasi-particle peak to higher energies diminishes the size of the D S Y loops 
involved in the coupled-channel problem. Therefore, the simultaneous dressing of the D and D s 
mesons in our self-consistent coupled-channel model produces a less bound A c (2595) resonance 
in nuclear matter, as already shown in Fig. [TJ From Fig. [2] we can see that the corresponding 
A C (2595)N~ 1 excitation mode of the D-meson spectral function appears approximately 40 MeV 
higher in energy than when only the D-meson dressing is considered. 

1=1/2, S=-l, C=-l 
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2860 2880 
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FIG. 3: Imaginary part of the I = 1/2, S = — 1 and C = — 1 D S N D S N scattering amplitude in nuclear 
matter at normal saturation density po an d zero temperature, as a function of the total energy Pq, for a 
total momentum P = and various approximations. 



The in-medium properties of the C = —1 mesons, D and D s , will be determined by the behavior 
of the corresponding DN and D S N amplitudes in the nuclear medium. In Fig. [3] we display 
the imaginary part of the D S N amplitude at normal nuclear matter saturation density and zero 
temperature as a function of the center-of-mass energy Pq, for various approximations: free (dotted- 
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FIG. 4: The same as Fig. [5] for the D (left panels) and D s (right panels) mesons. 

dashed line), Pauli blocking (dashed line), and the self-consistent calculation including only the 
dressing of D s (dotted line) and including both D s and D dressings (solid line). This is the most 
interesting of the two C = — 1 cases since the D S N system develops in free space a subthreshold 
bound state at 2906 MeV that couples significatively to D S N states. Therefore, this pole will be 
very sensitive to the medium effects. Indeed, when only Pauli blocking effects are considered, the 
pole moves about 40 MeV towards higher energy as expected. We observe very drastic changes 
when the dressing of the D and D s mesons is incorporated. The reason is that, as we will see, the 
in-medium quasiparticle peak of the D s meson experiences a strong attraction. This moves the 
in-medium threshold for D S N states below the position of the resonance, making its decay possible 
and quite probable due to the significant coupling to these states. 

The D and D s self-energies and spectral functions are shown in Fig. U] as functions of qo, 
including Pauli blocking effects, the additional self-consistent dressing of the given meson, and, in 
the case of the D s , incorporating also the dressing of the D meson in the DY intermediate states 
coupling to D S N . Again the thin solid lines indicate the q 2 . — m 2 ^ (left panel) and the q$ — m 2 ^ 
functions (right panel). 

The self-energy of the D mesons is quite smooth due to the absence of resonances in the DN 
sector. The repulsive character of the DN amplitude gives rise to a quasiparticle peak in the D 
spectral function appearing at higher energy than the D meson mass, by 30 MeV in the case of 
considering Pauli blocking effects only, or by 35 MeV when the additional self-consistent dressing 
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of the D meson is also taken into account. The repulsive mass shift obtained in Refs. 10ll . Il03l ] is 
somewhat smaller, of the order of 20 MeV for both cases. On the other hand, the self-energy of the 
D s meson shows a richer structure, which, in the case of Pauli-blocking, it is seen as a bump in the 
spectral function around 2000 MeV. This is a reflection of the presence, in the D S N amplitude, of 
a pole at 2906 MeV, which appears shifted about 40 MeV to higher energies when Pauli blocking 
effects are incorporated (Fig. [3]). The dressing of the D s meson smears this structure in such a way 
that one barely sees any trace of it in the corresponding spectral function. Moreover, the delta-like 
quasi-particle peak, appearing 60 MeV below the free D s mass when only Pauli blocking effects 
are considered, moves to slightly lower energies when the D s meson is dressed. Considering the 
additional dressing of the D meson in the related DY loops produces a substantial change in the 
D s self-energy. This is easy to understand from the results of Table HU where we see that the pole 
at 2906 MeV couples also very strongly to DY states. The loss of attraction in the region of the 
quasiparticle peak moves it towards a higher energy, and ends up being 50 MeV below the free 
mass and me rging with the resonant-hole strength. Our findings differ again quite strongly from 



those of Ref. 



101 



, which are dominated by an exotic coupled-channel molec ule a t 2780 MeV 



which is the equivalent to the pole at 2906 MeV found in the model of Ref. [105J and used in the 



lOlfl shows 



present work. As a consequence, the spectral function for the D s meson found in Ref. 
two distinct peaks, the quasi-particle one located about 10 MeV above the free D s mass, and a 
narrow resonance-hole mode located 150 MeV below. 
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FIG. 5: Imaginary part of the 1 = 0, DN — > DN scattering amplitudes in nuclear matter at po (left 
panel) and 2po (right panel), as functions of the total energy Pq, for a total momentum P = and two 
temperatures, T = (solid lines) and T = 100 MeV (dashed lines). The calculation in free space is also 
given for reference. 
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FIG. 6: The spectral function of the D meson at p = po (upper panels) and p — 2po (lower panels), two 
temperatures, T — (solid lines) and T = 100 (dashed lines) and two values of momentum, q — MeV/c 
(left panels) and q = 450 MeV/c (right panels). Dotted lines indicates the free mass of the meson 
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FIG. 7: The same as Fig. [5] for the D s meson 



In Fig. [5] we display the behavior of the A c (2595) resonance at two densities, po an d 2po, and 
two temperatures, T = and T = 100 MeV. We observe that the 35 MeV attraction felt by the 
resonance at p = po and T = MeV, gets substantially increased to almost 100 MeV at twice 
po - >■ 2po- The changes associated with temperature are also very significant. As already seen in 
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Ref. [103], increasing the temperature broadens the A c (2595) considerably due to the spreading of 
the D-meson strength because of the effect of Fermi motion. 

The effect of density and temperature in the spectral functions of the D, D s , D and D s mesons 
are shown in Figs. [6]to[9l where we compare results for two temperatures, T = (solid lines) and 
T = 100 (dashed lines), and two densities, p = po (upper panels) and p = 2po (lower panels), for two 
values of momentum, q = MeV/c (left panels) and q = 450 MeV/c (right panels), in the case of the 
complete self-consistent calculation, including Pauli blocking and dressing of mesons. A common 
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behavior in all spectral functions is that finite temperature moves the quasiparticle peak towards its 
free location. This is a reflection of the reduced size of the self-energy, because, being built up from 
an average over the smeared thermal Fermi distribution, involves higher momentum components for 
which the meson-nucleon interaction is weaker. Except for a few cases, increasing the temperature 
gives rise to wider quasiparticle peaks because of the increase of collisional widt h. H owever, the 
opposite effect is seen for the D meson in Fig. [6J As already discussed in Ref. [103] , this is due 
to the fact that the strength under this peak also receives contributions from £ C (2804)./V -1 hole 
excitations, which are washed out by temperature as any other resonant-hole mode. Consequently, 
the peak of the D-meson s pectr al function becomes narrower and more symmetric as temperature 
increases, similarly to Ref. 103 ]. 



The density effects observed in the spectral functions are also clearly understood. In general, we 
find that the self-energy roughly doubles its size when going from nuclear matter at normal nuclear 
matter saturation density to a system which is two times denser. This is consistent with the low 
density limit behavior and points at a weak density dependence of the in-medium meson-nucleon 
amplitude in this density region. This is the reason why, in general, the quasiparticle peak of the 
spectral functions at 2po are found approximately twice further away from the free space position 
and are twice wider than in the case of p$. 



IV. SUMMARY AND CONCLUSIONS 

We have studied the properties of open charm mesons, D, D, D s and D s , in nuclear matter at 
finite temperature within a self-consistent coupled-channel approach which uses, as meson-baryon 
interaction, a full i-dependent vector meson exchange driving force. 

The in-medium scattering amplitudes are obtained by solving the Lippmann-Schwinger equation 
at finite temperature including Pauli blocking effects, as well as D, D, D s and D s self-energies, 
paying a particular attention to their mutual influence. 

We have analyzed how our dynamically generated resonances are affected by density and tem- 
perature. As in other similar approaches, the resonances that couple strongly to intermediate 
states involving nucleons, move upwards in energy when Pauli blocking effects are considered, as 
a consequence of the loss of phase space. When the self-consistent dressing of the charm mesons is 
incorporated, the resonances gain attraction again. 

We have seen that dressing the D s meson has a non-negligible effect on the DN amplitude and 
on the properties of the D meson. Therefore, we conclude that a simultaneous in-medium treatment 
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of both mesons, as the one attempted in the present work, is necessary. Similarly, the in-medium 
properties of the D s and D mesons are interrelated and must be also considered together. 

The spectral functions of the D and D s mesons are quite rich. At T = MeV and normal 
nuclear matter density one finds a quasiparticle peak located below the corresponding free meson 
mass by about 50 MeV, as well as strength associated to resonant-hole excitations which, in the 
particular case of the D meson, is clearly visible as a narrow A C (2595)A^~ 1 excitation peak. 

In general, increasing the temperature has the effect of moving the quasiparticle peak towards 
its free location making it wider, as a consequence of a milder meson-baryon interaction and a 
larger amount of collisions. The exception found for the D-meson is naturally explained in terms 
of the mixing of the quasi-particle peak with a resonant-hole mode. 

For the densities explored, up to twice nuclear matter normal saturation density, we have found 
that the density effects follow the linear behavior expected for the low density regime: the self- 
energy roughly doubles its size when going from nuclear matter at normal saturation density to a 
system which is two times denser, indicating a mild density dependence of the in-medium meson- 
baryon interaction amplitudes. 

The enormous computational effort of the present work, which uses a coupled-channel formalism, 
an interaction that goes beyond the t — > limit, and the simultaneous consideration of the in- 
medium D and D s (D s and D) meson self- energies, has prevented us from incorporating the 
coupling to states involving vector-mesons. We are aware that, given the availability of models 
that permit dealing with these important degrees of freedom, our approach should be extended 
to the vector mesons such that it also includes, for instance, the D*N and D*Y channels in the 
C = 1, S = sector. We hope that, by first identifying which channels play a relevant role and 
which ones might be omitted, we can make progress toward this goal in the nearby future. 
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